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ABSTRACT 

Observations of high-redshift quasars show that the IGM must have been reionized 
at some redshift z > 5. If a source of radiation could be observed at the rest-frame Lya 
wavelength, at a sufficiently high redshift where some of the IGM in the line-of-sight 
was not yet reionized, the Gunn-Peterson trough should be present. Longward of the 
Lya wavelength, a damping wing should be observed caused by the neutral IGM whose 
absorption profile can be predicted. Measuring the shape of this damping wing would 
provide irrefutable evidence of the observation of the IGM before reionization, and a 
determination of the density of the neutral IGM. This measurement might be hindered 
by the possible presence of a dense absorption system associated with the source. 

Shortward of the Lya wavelength, absorption should be seen from the patchy 
structure of the IGM in the process of reionization, intersected in the line-of-sight. We 
show that a complete Gunn-Peterson trough is most likely to continue to be observed 
through the epoch where the IGM is partially ionized. The damping wings of the 
neutral patches should overlap if the proper pathlength through an ionized region 
is less than lh~^ Mpc; even in larger ionized regions, the characteristic background 
intensity should be low enough to yield a very high optical depth due to the residual 
neutral fraction, although occasionally some flux may be transmitted through large, 
underdense voids within an ionized region. The case of the He II reionization is also 
discussed, and we argue that helium was already doubly ionized by z = 3 throughout 
the IGM. 

The recently discovered afterglows of gamma-ray bursts might soon be observed at 
the very high redshifts required for these observations. Their featureless continuum 
spectrum and high luminosities make them ideal sources for studying absorption by 
the IGM. 



Subject headings: galaxies: formation - large-scale structure of universe - quasars: 
absorption lines 
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1. Introduction 

Observations of the spectra of high-redshift quasars blueward of the hydrogen Lya Une have 
demonstrated so far that the universe was already ionized by z ~ 5, the highest redshift at which 
any sources have been found so far (Schneider, Schmidt, k, Gunn 1991; Pranx et al. 1997). If the 
intergalactic medium (hereafter IGM) in the Unc-of-sight to a source was neutral, the resonance 
scattering at the wavelength of the Lya line should cause the "Gunn-Pctcrson trough" (Gunn & 
Peterson 1965), with an extremely high optical depth that would entirely absorb the flux from 
the source for any reasonable density of the IGM. At the same time, we expect that the known 
sources of ionizing photons (comprising quasars, stellar objects and hot gas) started to form much 
later than the recombination epoch, unless the amplitude of density fluctuations was much larger 
on small scales than predicted by all the models of structure formation. This implies the IGM 
should have been rcionized at some point. In fact, the COBE measurement of the electromagnetic 
spectrum of the microwave background (with a limit on the y parameter y < 1.5 x 10^^; Fixscn et 
al. 1996), as well as the presence of the Doppler peak in the intensity fluctuations, show that the 
IGM was reionized at z ^ 300 (Wright et al. 1994; Hu & Sugiyama 1994; Tegmark & Silk 1995). 

It is generally expected that the reionization was caused by the ionizing radiation from 
galaxies and active galactic nuclei. This implies that the IGM should have been reionized 
inhomogeneously, with every source ionizing first its immediate vicinity; the ionized bubbles would 
then fill a growing fraction of the volume in the universe until they overlap (Arons &: Wingert 
1972). Models of structure formation predict that reionization did not occur at extremely high 
redshift, but most likely at z ^ 30, given the epoch when the collapse of small-scale density 
fluctuations led to the formation of the flrst galaxies where enough stars could form to reionize 
the IGM (e.g., Tegmark, Silk, &; Blanchard 1994; Haiman, Rees, k, Loeb 1996). If a source could 
be seen at a redshift higher than that when most of the IGM was reionized, we should obviously 
see the Gunn-Peterson trough. And at redshifts were the reionization is already complete, we 
observe the Lya forest caused by the residual neutral fraction in the highly ionized IGM (in 
photoionization equilibrium with the background radiation), with density fluctuations originated 
in the gravitational collapse of gas in the developping large-scale structure. 

Thus, the question that naturally arises here is: how does the transition from one regime to 
the other take place? What should we expect to see in a spectrum where the increasing optical 
depth of the Lya forest with redshift turns into the Gunn-Peterson trough? This is the question 
that shall be addressed in this paper. Meiksin & Madau (1993) speculated that a new class of high 
column density absorbers might appear at redshifts higher than the epoch when the reionization 
was completed, due to the neutral patches left in the IGM. The possible observable signatures in 
the Lya spectra of sources seen at the time when the IGM still had this patchy ionized structure 
will be investigated in more detail here. 
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2. The Red Damping Wing 

We consider a source that is observed through a neutral IGM with uniform density. If the 
photons scattered in the Lja resonance were in a perfectly narrow line, we would expect to see 
a sudden drop in the flux at the Lya wavelength at the redshift of the source. In reality, the 
optical depth to Lya scattering is broadened both by the velocity distribution of the atoms and 
the natural width of the line. When the IGM is mostly neutral, the intrinsic line width is the 
dominant effect. The resulting absorption profile was discussed in Miralda-Escude & Rees (1997), 
but we reproduce this here for completeness. 

Any fluid element along a velocity interval dV along the line-of-sight (corresponding to 
the Hubble expansion velocity over a spatial interval dx = dV/ H{z)) has its optical depth 
spread in the observed spectrum according to r(F) = [rodV / c)Ra/'n: /[{V / cY + i?^]. Here, 
To is the Gunn-Peterson optical depth of the neutral IGM before broadening, and is given by 
TO = 2.1 X 10^ [l^b/i(l-F)/0.03] [Ho{l + zf''^/H{z)\ [{l + z)/%f/'^ , where Y is the helium abundance 
and H{z) is the Hubble constant at redshift z. The IGM is assumed to contain all the baryons, 
with density ilf, in units of the critical density. We have dcflned also Ra = A/(47rt'Q) = 2.02 x 10^^, 
where A = 6.25 x 10^ s~^ is the decay constant for the Lya resonance, and Va = 2.47 x 10^^ Hz is 
the frequency of the Lya line. As a result of this broadening, the observed absorption profile, at a 
wavelength separation AA on the red side of the Gunn-Peterson trough, should be given by 

= -r L/, (v/cY+Ri ■ 

Since tq ^ 1, the shape of this profile can only be observed when (A A) /A ^ Ra^ 

so we have 

r(AA) = Wtt ^— j =1.3x10 — (AA/A) . (2) 

Notice that the damping wing of the Gunn-Peterson effect has the optical depth falling as 
the inverse of the wavelength separation, as opposed to the inverse-square law for the damped 
absorption lines caused by high column density systems. The reason is that as we increase the 
wavelength separation, a larger pathlength through the neutral IGM contributes significantly to 
the optical depth, providing therefore a larger effective column density. 

If a source is found at very high redshift showing the Gunn-Peterson trough, the presence or 
absence of the damping wing with the predicted absorption profile can be an unambiguous test 
for the state of ionization of the intervening IGM. The observation of the Gunn-Peterson trough 
alone does not prove that a source is being observed behind neutral IGM, because the value of tq 
is so large that a small residual neutral fraction in a reionized IGM can reduce the transmitted 
fiux to undetectable levels. If the damping wing is observed, an obvious immediate application 
will be to infer the parameter tq from the observed width. At very high redshifts, we can use the 
approximation Q{z) ~ 1 independently of the value of ft at present, so H(z) ~ HqQ^/^{1 + z)^^^. 
Thus, by measuring tq we can infer the quantity QjjhCl"^^^. 
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There are two possible caveats in the measurement of this quantity. The first is that if 
the source is observed at the epoch when the IGM is being reionized, some regions along the 
line-of-sight should be mostly ionized and some others should be neutral. If the Hubble velocity 
across the HII regions is small compared to the width of the damping wing, the form of the 
damping profile should be essentially indistinguishable from the case of the homogeneous medium 
(and in that case, the width of the damping wing measures the space-averaged density of the 
fraction of the IGM that is neutral), but if the velocity is comparable then there will be random 
variations in the shape of the damping wing depending on the distribution of neutral gas close 
to the source. From equation ^ the width of the damping wing should be ~ 1500 Km s~^, or 
~ lh~^ Mpc, at z = 5. Notice that if there is only an ionized region around the source itself, with 
the rest of the gas in the line-of-sight that contributes significantly to the damping wing being all 
neutral, this does not affect the measurement of tq, but it only affects the redshift of the neutral 
gas closest to the source, which can be obtained from a two-parameter fit to the shape of the 
observed damping wing. 

The second possible problem is that the neutral IGM may not be homogeneous. In particular, 
any source of light is likely to be within a system that has gravitationally collapsed, and if the 
source has not ionized the surrounding gas producing an expanding cosmological H II region, we 
should expect it to be surrounded by a halo of accreting neutral gas. In that case, an additional 
column density of neutral hydrogen is added, very close to the redshift of the source. One could 
always try to circumvent the problem by adding the "local" column density as a third parameter 
to the fit to the damping wing profile, but for a realistic accuracy of the spectrum the errors will 
of course greatly increase as more parameters are added to the fit. 

As an example, we plot in Figure 1 the absorption profile on a source at z = 9 {solid line), 
for To = 4.3 X 10^. We have used the exact expression for the profile of the damping wing due to 
the neutral IGM, presented in the Appendix, where we do not assume (A A) /A ^ 1 (as in eq. |2|), 
and we include the effect of the cosmological evolution of the Gunn-Peterson optical depth, and 
the exact expression for the cross section of the Lya resonance line. The correction on the optical 
depth relative to the approximation in eq. § is about 20% when AA/A = 0.02, and reaches a factor 
of 2 at AA/A = 0.06; it will therefore be important to use the exact expression in the Appendix 
for practical applications. We also plot in Fig. 1 the result of adding to the IGM optical depth an 
absorbing system with column density Nhi = 2 x lO^'^ cm~^ at the source redshift z = 9 (dotted 
line) . Finally, the dashed line shows again the IGM optical depth without any additional absorbing 
system, but with the parameters z = 9.009, tq = 4.9 x 10^, which provide the best fit to the dotted 
line. It is clear from this example that even with an accuracy of 1% in the determination of the 
shape of the damping wing, there is still a large error in the determination of tq if an absorber of 
unknown column density is present at the redshift of the source. However, one may be able to 
determine the source redshift independently (for example, from associated metal-line absorption), 
in which case the ambiguity between the hydrogen column density of an associated absorber and 
To should be much reduced. 
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3. Flux Transmission on the Blue Side of the Gunn-Peterson Trough 

Moving down in redshift from the source, a line-of-sight will characteristically traverse several 
cosmological H II regions and neutral patches over the redshift interval corresponding to the epoch 
when the IGM was being reionized. The process of reionization can be very complicated: the 
sources of ionizing photons may have a wide range of luminosities and may be short-lived, so some 
ionized regions will start recombining if their source fades, becoming partially neutral until they 
are ionized again by another source. The sources could be highly clustered, so the ionized regions 
may start being very small and due to faint, individual sources, and later grow to a much larger 
scale and be ionized by clusters of sources, while large neutral patches may still remain in the IGM. 
And of course, the HII regions expand into a highly inhomogeneous IGM which is gravitationally 
collapsing into new large-scale structures. The case of a source of constant luminosity surrounded 
by a homogeneous medium can be solved analytically (Shapiro & Giroux 1987). 

An ionized region intercepted by the line-of-sight over a total length Li = Vi/H{z) may cause 
a "gap" in the Gunn-Peterson trough if the damping wings of two neutral patches on each side do 
not completely overlap. Taking the parameters used in Fig. 1, we see that if we require that at 
least a fraction = 0.36 of the flux be transmitted in the middle of the gap, then the length of 
the intersected ionized region must be at least Vi/c > 10~^ (because at a separation V/c = 5x 10~^ 
from the edge of a neutral zone, the optical depth of the damping wing in Fig. 1 is r = 0.5). 
Since tq oc (1 -|- z)^^"^ , this minimum size scales with redshift as Vi/c > 10^^ [(1 -|- z)/10]^/^ in the 
approximation of eq. |2|, or a proper length Li > lh~^ Mpc (these scales are also proportional to 
0,bh{l — Y), but we are setting this quantity to 0.03 in this section). 

In order to observe any transmitted flux to the blue of the Gunn-Peterson trough on a source 
at redshift Zs, some HII region of this size has to be present on the line-of-sight at a redshift 
1 + Zi > 1 + zp = {1 + Zs)/l.l8, since at lower redshifts the Ly/3 Gunn-Peterson trough blocks 
the flux anyway. If the last neutral patch in the line-of-sight is located at Zn > zp, then flux may 
be observable between the blue Lya damping wing due to this last neutral patch and the Ly/3 
damping wing of the neutral IGM close to the source, and some additional small gaps may be 
observable due to large HII regions at z > Zn- li Zn < zp, then only these gaps due to isolated HII 
regions may be observable. 

The existence of ionized regions with a proper size larger than lh~^ Mpc for a substantial 
period of time before reionization was over depends, of course, on the type of sources causing the 
ionization. Large HII regions would be expected if the radiation was dominated by very luminous 
sources; the luminosity required to reach a proper size Mpc at z = 10 is similar to that of the 
most luminous quasars known. Reionization by early galaxies has been considered more likely in 
the usual hierarchical models for structure formation (Tegmark, Silk, & Blanchard 1994, Haiman 
& Loeb 1997, MR97). In this case, large HII regions might still exist if the sources were highly 
clustered, probably in the sites of formation of the future galaxy clusters. Early galaxies might in 
fact be clustered on large scales due to a high biasing factor, which is expected if galaxies form in 
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density peaks and when the power spectrum has a slope close to n = —3, as is generally predicted 
to be the case on small scales (Bardeen et al. 1986). On the other hand, if reionization is caused 
by low-luminosity sources which are not highly clustered on a proper scale of lh~^ Mpc, then by 
the time the ionized regions had overlapped to this scale there should be very few neutral patches 
remaining in the IGM, (except for the neutral regions that are dense and self-shielded, i.e., regions 
which are neutral because of their high recombination rate, rather than the lack of sources in their 
vicinity) . 

Whether any flux is actually seen on the blue side of the damping wing of the last neutral 
patch in a line-of-sight, or in any gap due to an isolated ionized region, depends of course on 
the neutral fraction in the ionized medium being low enough. The neutral fraction should be in 
ionization equilibrium with the radiation field from the sources in the H II region, and will cause 
a "Lya forest" , as observed up to z ^ 5 so far; even a very small neutral fraction can completely 
absorb any remaining flux that is left between the damping wings of the neutral patches. 

To calculate the absorption in the ionized regions, let us parameterize the mean volume 
emissivity of ionizing photons as 

6 = ^iV. , (3) 

where Ue is the mean electron density in the ionized IGM and t is the age of the universe. Thus, 
when Nr = 1 the emissivity is such that one ionizing photon will be emitted for each electron 
over the age of the universe. Obviously, at the end of reionization the factor Nr should be greater 
than unity, both because the comoving emissivity should increase with time (and therefore the 
integrated number of emitted photons is less than et), and because electrons can recombine several 
times during reionization, increasing the number of photons that need to be emitted. The photon 
intensity is typically e Li/ (Sir) (where we take the typical mean-free-path of photons to be ~ Li/2), 
and the implied photoionization rate is F = eLia/2, where a is the photoionization cross section 
averaged over all emitted photons; for a typical power-law slope Fy oc v"^'^ near the hydrogen 
ionization edge, ~ 2 x 10~^^ cm^. In an HII region having the minimum size that allows flux to 
be transmitted (using as before the condition that the optical depth between the two damping 
wings reaches a value less than unity), the length is Li = [c/ H{z)]{'iTQRa/'K) (see eq. The 
neutral fraction in a region of average density is equal to aUe/T, and the optical depth in the 
ionized region due to this neutral fraction is: 

= -TV, - = = 4 

(Ue/t) Nr (ToRa/n) [c/H{z)] <7 3caRa Nr Nr ' ^ ' 

Here, T4 is the temperature of the ionized gas in units of 10^ K; the temperature dependence is 
due to the recombination coefficient. Notice that Ti is the optical depth in the HII region if the 
gas density is equal to the mean baryonic density. In reality, the gas density should of course 
fluctuate and the optical depth should consequently vary with wavelength over the interval of Lya 
absorption of the HII region. 

Thus, we have found the result that the optical depth to Lya scattering in an ionized region 
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that is just large enough to emerge from the damping wings of the Gunn- Peterson trough is 
independent of redshift and of any cosmological parameters, and depends only on the dimensionless 



Let us explore more carefully the interpretation of A^^. There are two different cases that 
need to be considered. In case I, the sources of ionizing photons turn on quickly enough so 
that recombinations are negligible compared to the number of photons that need to be emitted. 
Therefore, the IGM is reionized when one ionizing photon has been emitted for each electron, 
and Nr = e/e, where e is the time-averaged comoving emissivity. Unless the sources are highly 
synchronized, Nr will not be larger than ~ a few. The other possibility, case II, is that reionization 
is limited by the rate at which electrons can recombine. In that case, Nr is simply the mean 
number of recombinations for each electron that take place over the age of the universe, t: if 
the sources turn on over a time tr, each electron undergoes Nr{tr/t) recombinations during the 
time tr, and the required emissivity is e = {ne/tr)Nr{tr/t) = {ne/t)Nr. This mean number of 
recombinations is 



where the last term is the electron clumping factor. Thus, at a fixed redshift, the only way for 
Nr to be large is that the IGM be highly clumped, and the clumps need to have a large covering 
factor through a line-of-sight in an ionized region in order to absorb the photons effectively (if 
the covering factor is small, the condition that the mean recombination rate is equal to the mean 
photon emission rate implies that the clumps are mostly neutral, therefore not increasing the 
electron clumping factor). But if gas clumps with a high covering factor are responsible for the 
recombinations and, therefore, for the absorption of most photons before they reach the edge of the 
H II region and contribute to increase its size, then our estimate of the photoionization rate in the 
HII region is not valid because we assumed the mean- free-path of the ionizing photons to be Li/ 2, 
and it should now be reduced by the covering factor of the clumps. As long as the recombinations 
do not take place mostly at the edges of the HII regions (due to a higher gas density there), 
Nr cannot be larger than a few even for a very clumpy IGM. Notice that A',, is very large when 
z ^ 10, but of course the required size of an HII region Lj to allow for flux transmission between 
the two damping wings becomes more implausibly large as the redshift increases. 

We therefore conclude that the optical depth Tj of an ionized region before reionization is 
complete should generally not be smaller than ~ 100. Rare exceptions to this are nevertheless 
allowed, because the photoionization rate inside the HII regions that we estimate above is only 
a mean value; for example, a line-of-sight may pass very close to a very luminous source. At 
the same time, the optical depth Tj is the value expected when the IGM density is equal to the 
mean baryonic density of the universe, but in reality a " Lya forest" should be caused by the 
IGM density fluctuations, just as observed at lower redshift. In photoionization equilibrium, the 
optical depth observed at a given wavelength is proportional to the square of the gas density. 
Therefore, the optical depth should be reduced by a factor of 100 in a void that is underdense 
by a factor of 10, allowing for some flux to be transmitted. The transmission of any flux when 



number A',., defined in eq. |3|. 




(5) 
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Ti is so large requires a void that is not only sufficiently underdense, but also sufficiently large 
to prevent the thermal broadening of the absorption due to the gas in adjacent structures, with 
typical overdensities near unity and temperature ^ 10"^ K, from overlapping and obstructing the 
window. For example, for absorbers with central optical depth ~ 100 and velocity dispersion 
(7 = 10 Km s~^ the thermal wings should extend to about Sa on each side, requiring the void to 
have a diameter of at least ^ (60Kms'~^)/iJ(2;). 

Possibly, the epoch of partial ionization in the IGM may be better observed in the LyP 
line (and higher order Lyman series), where the optical depth is reduced everywhere by a factor 
0.16 relative to Lya . This could work in a source observed at a redshift close to the end of the 
reionization epoch, and if the intensity of the ionizing background rises quickly enough that, by 
the time the scale factor of the universe has increased by 1.18, significant flux is already being 
transmitted through the Lya forest. But of course, contamination by the Lya forest will make the 
analysis of the Ly/5 spectrum very complex. 

To summarize, the picture that emerges from this discussion is the following: at the redshift 
Zn of the last neutral patch of the IGM on the line-of-sight from the source, a blue damping wing 
of absorption should be produced by the neutral gas. However, the ionized IGM at redshifts 
immediately below Zn should normally have a large enough neutral fraction, in equilibrium with 
the ambient intensity of ionizing photons, to completely obstruct the flux from the source at the 
Lya resonance, thus preventing the observation of the shape of the blue damping wing. The 
highest redshift at which transmitted flux will be seen should be that of a sufficiently underdense 
and large void, which can open up a "window" in the saturated absorption of the ionized IGM. 
Some of these voids might be located in the region where the optical depth of the blue damping 
wing is still significant, or even between two damping wings encircling an ionized region at 2; !> Zji 
with a size larger than Lj. 

At z < Zn, the intensity of the ionizing background, Jhi, should rise as the mean-free-path 
of the photons increases, leading to an increase in the number of "windows" of transmitted flux 
associated with voids. The redshift evolution of these "windows" will provide information on 
the evolution of Jhi- How fast should the rise of Jhi be? This depends on the rate at which 
the ionized material can recombine. In case I (where the mean rate of recombinations is much 
lower than the mean rate of emissions before the HII regions overlap), reionization occurs after 
one ionizing photon has been emitted for every proton in the IGM. At this point, the photon 
mean-free-path should rapidly increase by a large factor starting from the typical size of the 
HII regions, thereby causing a fast increase in the fraction of transmitted flux in Lya spectra. 
The intensity of the background rises proportionally to t — ti, where ti is the time at the end 
of reionization. The increased background radiation penetrates regions of denser gas, increasing 
the global rate of recombinations in the universe until they can balance the emission rate. This 
stops the linear growth of Jhi with time, and the photon mean-free-path has grown to the mean 
separation along a random line-of-sight between the Lyman limit systems that have now been 
ionized, where enough recombinations take place to match the mean emission rate. The subsequent 
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evolution is known observationally: at z ~ 4 to 5 the abundance of Ljnnan limit systems is still 
large enough to be the main factor limiting the intensity Jhi, and by z ^ 2 the mean- free-path 
has grown to the horizon scale, and Jhi is limited mainly by the redshifting of photons. 

On the other hand, in case II recombinations are already important before the H II regions 
start to overlap. At this time, the mean-frcc-path between Lyman limit systems associated with 
dense, self-shielded structures is about the same as the typical size of the HII regions. Thus, there 
is no sudden increase of the mean-free-path when the HII regions overlap. As the emission rate of 
photons increases and the universe expands, the self-shielded regions shrink in size and the newly 
ionized dense gas enables the global recombination rate to keep up with the emission rate, and 
this continues into the observed epoch at z < 5 as discussed previously. This leads to a much 
more gradual increase of the mean-free-path and of J^j compared to case I. Thus, the width of 
the redshift interval over which the Lya forest increases its opacity, leading to the Gunn-Peterson 
trough, should give us information on the dumpiness in the IGM, which determines the rate of 
recombination at various stages of the reionization. 



4. The Hell Gunn-Peterson Trough 

The first observations of absorption by Hell at A = 304A have been performed over the 
last few years (Jakobsen et al. 1994; Tytler et al. 1995; Davidsen, Kriss, & Zheng 1996; Hogan, 
Anderson, & Rugers 1997; Reimers et al. 1997). The observation by Jakobsen ct al. (1994) was 
consistent with a complete "Gunn-Peterson trough", although the signal-to-noisc that could be 
achieved allowed them only to place an upper limit to the mean transmitted flux of 20% of the 
total at z ~ 3.2. The Hell Lya optical depth from a uniform IGM where all the helium is only 
once ionized is 

ro,Heii = 2.3 X 10^ {QhhY/OM) [Ho{l + zf/^Hiz)] [(1 + z)/Af^ . (6) 

Therefore, the observation of Jakobsen et al. could imply that the helium in the IGM was not 
yet doubly ionized, but it was also consistent with a typical Hell fraction much smaller than 
unity (this remains true in a realistic model of an inhomogeneous IGM with on-going structure 
formation). 

Davidsen et al. (1996) were the first to find evidence for the presence of transmitted flux to 
the blue of the Hell Lya line in another quasar, with a fraction 0.36 =b 0.03 over the redshift 
range 2.2 < z < 2.7. The detection of flux immediately proves that the helium in the IGM was 
mostly twice ionized at z < 2.7, although it is possible that some patches of Hell still remained 
at this epoch if the double reionization of helium was not yet complete. In fact, Reimers et al. 
(1997, hereafter R97) detected transmitted flux with much higher resolution in a second quasar, 
QSO HE 2347-4342, in the range 2.8 < z < 2.9. All the transmitted flux appears in two windows, 
one at A = II6OA of width 4rmA, and another at A = 1174A of width 2A. The corresponding 
HI Lya absorption spectrum (see Fig. 4 in R97) shows some weak absorption features in the 
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wider window, and no detectable absorption in the narrower one. In particular, H I absorption at 
A = 1159 A, in the middle of the wider Hell window, is clearly detected. On the other hand, there 
are other regions free of any detectable absorption in HI (e.g., at A = 1171 A) where the flux is 
completely absorbed in He II . This demonstrates that the various regions in the IGM observed in 
this line-of-sight cannot all be photoionized by an ionizing background with approximately the 
same intensity, as shown in detail by R97 (see their Fig. 5). It is clear that, while in the two 
windows of the He II spectrum the helium in the IGM must be mostly twice ionized, the other 
regions with low H I absorption must either have most of their helium only once ionized, or must 
have a much lower intensity of Hell ionizing photons. 

The double ionization of helium could very well occur at a later epoch than the hydrogen 
ionization. If recombinations are not important during the reionization (case I as defined in 
the previous section). Hell will be ionized later if the ratio lueii /Ih i of emitted photons 
above the Hell ionization edge and above the HI ionization edge is lower than 0.08. But if 
many recombinations take place (case II), then we only need IhsII /Ihi < 0.5, because Hell 
recombines 5.5 times faster than HI (Miralda-Escude k, Rees 1994). If quasars or starburst 
galaxies produce the ionizing background, then IhcII /Ih I ^ 0.15 (Haardt & Madau 1996 and 
references therein). Since even in the absence of dumpiness, the Hell recombination time is equal 
to the Hubble time at z ~ 3, a delayed reionization of Hell relative to HI is clearly not just 
possible, but very likely. 

This suggests that we could be observing the kind of patchy helium medium that we discussed 
previously for the case of hydrogen, and that the two windows of transmitted flux observed by R97 
could be He III regions surrounded by Hell patches. However, the story for the Hell ionization 
is quite different from the case of hydrogen. The difference lies in the much larger thickness of 
the ionization fronts in the case of helium, due to the lower helium abundance by a factor 0.08, 
and the lower cross section of Hell by a factor 4. At z = 3 and in a region of density equal 
to the mean, a He HI ionization front should have a proper thickness ~ 1 Mpc, comparable to 
the largest He III regions that could be ionized by a single quasar. The corresponding velocity 
width is ~ 500Kms~^ (much larger than the expected width of the damping wing for the Hell 
Gunn-Peterson trough, which is smaller than that of HI by a factor 0.08 at the same redshift). 

The large thickness of the Hell ionization fronts implies that, during the epoch when the 
low-density IGM still contains a high fraction of He II , it is not likely that any windows of 
transmitted flux can be seen through the Hell Gunn-Peterson absorption. In fact, if we repeat 
the analysis of the previous section to obtain the optical depth Ti^jje II of a He HI region due to 
the Hell remaining in photoionization equilibrium, but this time expressing the answer in terms 
of the length of the He III region Li = Vi/H{z), we obtain 

To,Hen 750Kms~^ 
Ti,Hen=^^ Vi • 

Since the windows of transmitted flux observed in R97 have velocity widths similar to 

Vi = 750Kms~^, and Nj. cannot be a large number due to the same arguments used in the last 
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section, we conclude that the optical depth in a He III region would in fact be close to tq^hs ii (eq. 
^) if the double reionization of helium was not yet complete. Another simple way to understand 
this result is that when the thickness of the ionization front is similar to the size of the He III 
regions, the fraction of Hell in the He HI regions is not much less than unity. 

The fact that the Davidsen et al. observations can be fitted with models where the helium 
is already doubly ionized ai z = 2.7, with the expected shape of the spectrum for the ionizing 
background when quasars are the dominant sources (Miralda-Escude et al. 1996, Croft et al. 1997) 
also suggests that the Hell reionization should be complete, and the state of the IGM should 
not be radically different at z = 2.8. The obvious differences in the Hell to HI column densities 
presented by R97 can be explained by large fluctuations in the intensity of the background above 
the Hell ionization edge, which are expected due to the large number of optically thick absorbers 
in He II , but they do not demand a large He II fraction in the regions free of H I absorption with 
no transmitted flux at Hell . As an example, consider a typical void in the IGM at z ~ 3 with a 
gas density a few times below average, where the hydrogen neutral fraction is ~ 10^^. This yields 
an optical depth in the HI Lya spectrum of ~ 2% (e.g.. Croft et al. 1997), which would not be 
detected in any of the regions that are apparently free of absorption in the H I spectrum in Fig. 4 
of R97 due to the uncertain continuum level. Given the mean spectrum of the ionizing ackground 
obtained from quasar models (Haardt & Madau 1996), with Jh i /Jiieii — 50, the fraction of 
Hell in the same void should be 10~^, yielding a Hell Lya optical depth of 0.4, which could 
perfectly well be the case for the two windows of flux in R97. But in another similar void where 
Jh I / J He II = 500 (due to absorption of the nearest quasars), the Hell optical depth should be 4 
while the Hell fraction is still as low as 10~^. 

5. Conclusions 

In this paper we have attempted to predict the observable features that should be present in 
the spectrum of a high-redshift source seen behind a region of neutral IGM, still to be reionized. 

The first feature is the damping wing of the Gunn-Peterson trough on the red side of the 
Lya line. The width of the absorption profile of this damping wing provides a method to measure 
the density of the neutral IGM near the redshift of the source. There are two main difficulties 
for this measurement: the possibility of an absorption system large enough to contaminate the 
profile of the damping wing, and uncertainties in the intrinsic emission spectrum of the source. 
The second difficulty is particularly severe if the emitting source is a quasar, since quasars have 
strong, broad Lya emission lines with profiles that are highly variable. The second feature we 
have discussed is the transition from the Gunn-Peterson trough to the Lya forest. We have found 
that the Lya forest should gradually become thicker with increasing redshift as the epoch of 
reionization is approached; at redshifts where the IGM was only partially ionized, the presence of 
any transmitted fiux in individual cosmological HII regions surrounded by neutral patches of the 
IGM should be rare, due to the damping wings of the neutral zones and the very high Lya optical 
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depth that is still caused by H II regions with the characteristic intensity of the ambient ionizing 

background that should be expected. It is also possible that some transmitted flux could be seen 
from HII regions in Ly/3, although that would of course be contaminated by the Lya forest. 

What type of sources can we expect to discover at ever higher redshifts in the future, to 
finally be able to probe the epoch of rcionization? Quasars are the most luminous known sources, 
but there are signs that their abundance is declining at z ^ 3 (Warren, Hewett, Sz Osmer 1994; 
Schmidt, Schneider, & Gunn 1995) and it is not clear if they can be found at much higher redshifts. 
Young galaxies must necessarily have existed when the IGM was partially ionized if they are the 
sources that caused the reionization, but the expectation is that they are hopelessly faint for 
obtaining good quality spectra (Haiman & Loeb 1997; MR97). Recently, a new class of luminous 
extragalactic sources in the ultraviolet has been identified: the afterglows of gamma-ray bursts 
(Costa et al. 1997a; van Paradijs et al. 1997). The afterglow was predicted from the cosmological 
fireball model (Meszaros Sz Rees 1997; Vietri 1997), and found to be in good agreement with 
observations (e.g., Wijers, Rees, & Meszaros 1997a; Waxman 1997). The third observation of an 
afterglow (Costa et al. 1997b; Djorgovski et al. 1997) resulted in the detection of an absorption 
system in the spectrum at z = 0.83 (Metzger et al. 1997), which has definitely put to rest any 
doubts on the cosmological nature of GRBs. The optical magnitudes of the GRB afterglows are 
similar to the luminous high-redshift quasars. It is quite plausible that these afterglows will be 
found in the coming years at redshifts higher than any other known sources; in fact, if GRBs 
are produced in objects belonging to young stellar populations, then they should take place soon 
after the first stars were formed in the universe (Wijers ct al. 1997b). At very high redshifts, the 
probability of gravitational lensing is high, and this can help in boosting the apparent brightness 
of some GRBs and their afterglows. The very important advantage of GRB afterglows is that their 
intrinsic spectrum is supposed to be a featureless power-law on small wavelength ranges, because 
the radiation originates from synchrotron emission in a relativistic shock. This eliminates the 
severe systematic error of measuring the absorption profile of a damping wing near an emission 
line in the case of a quasar. Thus, unless gamma-ray bursts did not take place for some reason 
until a much later time than the epoch when the first stars (or the first sources that caused the 
reionization) formed, their afterglows should become a precious tool to study the IGM at very 
high redshift. 

I am grateful to Martin Rees for numerous discussions that resulted in many of the ideas 
presented here. I also thank Peter Meszaros and Eli Waxman for discussions on gamma-ray bursts. 
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Appendix 

We calculate here the accurate profile of the damping wing of the Gunn-Peterson trough 
caused by a homogeneous neutral IGM. The scattering cross section of the Lja resonance line by 
neutral hydrogen is given by (see Peebles 1993, §23) 



^ ' Svr {oo - coaf + (A2/4) (cv/uJaf ' ^ ' 

The notation is the same as that used in §2, and LOa = ^irva = 2Trc/Xa- The second term in 
the denominator can be neglected, because we are only interested in the shape of the damped 
profile far from the center of the line (i.e., when ||a; — (Ja\\ » A, where the resulting optical depth 
is not very large. The term in the numerator causes the classical Rayleigh scattering, and was 
neglected in §2 assuming — Wa|| coa, but this is a poor approximation since the damping 
wing can be quite broad, and wc keep this term here. We assume that the IGM has a constant 
neutral hydrogen comoving density uq at redshifts Zn < z < Zg, where Zs is the redshift of the 
source, and the density is negligibly small at z < Zn. The optical depth observed at a wavelength 
A = Aa(l + Zs) + AA is 

where 5 = AA/[Aa(l + ^s)]. The Gunn- Peterson optical depth outside the damping wing is 
To{zs) = 3A^Ano/[87rciI(zs)], and using H{z) oc (1 + z)^/^, we obtain: 



tqRc r» dz fl + z\^^/'^ , ^^_4 



l + z 



{l + Zs){l + 5) 



(3) 



This reduces to: 



r(^X) = '-^il + Sfl- r^,, (4) 
TT Jxi {l-xy 

where xi = (1 + Zn)/[(1 + -Za)(1 + S)], X2 = (1 + and the result of the integral is 
[ dxx^/^ a;9/2 9 7/2 9 5/0 , 3/0 ^1/2 9, 1 + x^/^ 
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Fig. 1. — Absorption profile of the damping wing of the Gunn-Peterson trough, predicted using 
the equation given in the Appendix. The fraction of transmitted flux is shown as a function of the 
wavelength interval AA from the Lya wavelength of the source, Aq(1 + Zs). The solid line assumes 
the parameters Zs = 9, tq = 4.3 x 10^, and z„ = 7 (where the IGM is assumed to have a constant 
neutral density at Zn < z < Zg, and be fully ionized at z < Zn)- The dotted line is the case where 
an absorption system with Nhi = 2 x lO^*^ cm~^ is present at Zg = 9, which could be due to a 
galaxy where the source is located. The dashed line shows the damping wing profile without any 
additional absorption system, but with the parameters changed to Zs = 9.009 (although AA is still 
defined as the wavelength interval from the Lya line at z = 9), tq = 4.9 x 10"'', and provides the best 
fit to the dotted line assuming there is no absorption system at the source redshift. This illustrates 
the difficulty in measuring tq in the presence of an absorption system. 



